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O ! Abstract 

Q.^ I We put bounds on the variation of the value of the fine struc- 

■^ ' ture constant a, at the time of Big Bang nucleosynthesis. We study 

p . carefully all light elements up to ^Li. We correct a previous upper 

' I limit on |Aa/a| estimated from ^He primordial abundance and we 

^H ' find interesting new potential limits (depending on the value of the 

c/3 . baryon-to-photon ratio) from ^Li, whose production is governed to a 

large extent by Coulomb barriers. The presently unclear observational 



S 



situation concerning the primordial abundances preclude a better limit 
^ ■ than |AQ/a| < 2 • 10"'^, two orders of magnitude less restrictive than 

previous bounds. In fact, each of the (mutually exclusive) scenarios 
of standard Big Bang nucleosynthesis proposed, one based on a high 
value of the measured deuterium primordial abundance and one based 
on a low value, may describe some aspects of data better if a change 
in a of this magnitude is assumed. 

1 Introduction 

Physicists have long speculated (at least since the time of P. A. M. Dirac 
[jl]]) about possible variations of the fundamental physical constants. The 

*E-mail addresses: lbe@physto.se, sergio@physto.se, rub@physto.se. 



fine structure constant, a = e'^/in, is especially interesting to test, being 
dimensionless and accurately known experimentally.^ 

Several attempts to constraint the time variation of a have been made 
in the last years [0, |^, ^, ||. The methods involved in these computations 
are quite different and the results are complementary since these calculations 
limit the variations of a at different cosmological times. In the absence of 
a particular model for the time dependence of a, there is no compelling 
reason to fit it in any particular way. Thus, it is important to find limits on 
the variations of a at different epochs. In particular, there may be a strong 
dependence on the cosmological redshift parameter z which could make most 
direct measurement methods insensitive, whereas an indirect method like the 
one we are investigating - primordial nucleosynthesis - puts us nearer the 
epoch of the unknown physics of the Big Bang itself where certainly the 
Standard Model of particle physics is inapplicable. 

Despite the lack of explicit models for the time variation of a, it should be 
kept in mind that there exist arguments, e.g., from fundamental theories with 
extra dimensions [0) H? which support the general idea of non-constancy 
of the low-energy parameters. However, models in which the coupling is 
governed by some condensate that varies in time are not yet fully understood. 

In some (e.g., p|) the dynamics has a "pre Big Bang" phase. After infla- 
tion, the universe is such that the dilaton is already massive and the couplings 
are set at the bottom of the potential without much further variation. 

In another class of models the dilaton remain massless but it decou- 
ples from matter. However, T. Damour and A. M. Polyakov predict possible 
effects of dilaton-induced changes in a and other quantities relevant to pri- 
mordial nucleosynthesis to be of the order 10^®, though in this model other 
dilaton effects are more important. In view of the primitive state of theory for 
an eventual time variation of fundamental couplings, we use a phenomeno- 
logical approach and assume a different value of a than the present one at 
the time of nucleosynthesis and investigate what the observable consequences 
would be. 

Direct measurements in the laboratory have given a limit on the variation 
|Aa/a| < 10~^^ over a period of 140 days [Q. Astrophysical observations of 
spectra of high red-shift quasar absorption lines have given limits of | Aa/aj < 
10~^ — 10~^ for z ranging from 0.25 to 3 0. The geological limit from 

^The presently established value of a at zero momentum transfer is given by 
a~i = 137.035989(61) [|. 



the Oklo natural nuclear reactor is about |Aa/a| < 10^^ over a period of 
1.8 billion years 0. Recently, it was argued that from future observations 
of fluctuations in the cosmic microwave background radiation, variations of 
a could be bound by |Aq;/q;| < 10~^ — 10"^ for z ~ 10^ |jlO|. Finally, 



assuming a particular model for the a-dependence of the neutron-proton 
mass difference, a limit can be extracted from the ^He primordial abundance 
for z ~ 10^ - lO^o [|. 

This work deals with contraints from nucleosynthesis, which have the 
advantage of probing the earliest cosmological epoch where data exist and 
where the basic physical processes are known from laboratory experiments 
work. Our analysis corrects that of and extends it since we have not only 
considered the '^He abundance but also the abundances of other light nuclei. 

Besides consituting a more thorough analysis of the nucleosynthesis bounds, 
our treatment has the advantage of bypassing the difficult theoretical prob- 
lem of the dependence of the neutron-proton mass difference on a, relying 
on the fact that the other abundances are only weakly sensitive to this mass 
difference. The abundance of the other light elements as a function of a can 
be extracted with smaller ambiguity, and the dependence is in addition very 
steep, in particular for ^Li (at least in certain ranges of the baryon-to-photon 
ratio, r]). Of course, at the present time these advantages are to a large extent 
balanced by the disadvantage of observationally less well known abundances, 
something that may hopefully improve in the future. 

In the next Section we review the most important aspects of the Standard 
Big Bang Nucleosynthesis (SBBN) model. Then, in Section 3 we explain how 
a enters into the SBBN scenario. In Section 4, we study the a-dependence of 
the relevant nuclear reactions involved in primordial nucleosynthesis. Finally, 
in Section 5 we estimate the corresponding limit placed on |Aa/a| at 2; ~ 
10^ — 10^°. We find that, given the present observational uncertainties of 
light-element abundances, it is not possible to put a more stringent bound 

than 

lAal o / , 

^< 2- 10-2, 1) 

a 

two orders of magnitude less restrictive than the one claimed in |^. 



2 The fine structure constant in the SBBN 
scenario 

In order to give a complete account of the effects of a variation of the fine 
structure constant on the primordial abundances, we will briefly review the 
SBBN model. We will focus our attention on those aspects concerning the 
evolution of the abundances of light nuclei. 

The production of light elements according to the SBBN model is com- 
monly divided into three stages [111], O- 



• First Stage: Statistical equilibrium (T ^ 1 MeV; t ^ 1 sec) 

During this first stage (as in the other two), the universe is radiation- 
dominated. The relativistic degrees of freedom are photons, electrons, positrons 
and the three light neutrino species. The weak interactions that interconvert 
neutrons and protons, i.e., 

n ^ p + e~ +Ve (2) 

n + Ue^ p + e~ (3) 

and 

n + e^ ^ p + T'e (4) 

are rapid enough to keep them in statistical equilibrium. The neutron-to- 
proton ratio is, then, given by its equilibrium value (we use units such that 
the speed of light, Planck's constant and Boltzmann's constant are all set to 
unity) : 

(I) =''^^'^ (') 

where Am = uin — rup is the neutron-proton mass difference. 

Moreover, at this temperature, not only are the rates for the weak reac- 
tions (H)-(§) faster than the universe expansion rate, but so are the nuclear 
reaction rates responsible for the production of light elements. Light nu- 
clei are then both in kinetic and chemical equilibrium or nuclear statistical 
equilibrium (NSE), their corresponding abundances being given by |TT[] 

Ya= ^A[C(3)^-V(i-^)/22(3^-5)/2]Ai(^)'^^-'^/' X 

V^-'YfY^^-^ exp{BA/T) (6) 



where t] is the present baryon-to-photon ratio, rriN is the nucleon mass, Ba 
is the binding energy, qa counts the number of degrees of freedom of the 
nuclear species {A, Z) and ( is the Riemann zeta function.0 

Due to the fact that the number density of photons is so large relative 
to that of baryons {rj ~ 10~^°), the abundances of composite elements are 
completely negligible and their synthesis does not truly start at this epoch. 

• Second Stage: Neutron-proton "freeze-out" (T ~ 0.8 MeV; t ~ 2 sec) 

After neutrino decoupling, at about the time electron-positron pairs anni- 
hilate, the second stage of primordial nucleosynthesis takes place. The weak 
reactions (^-(^ become slower than the expansion rate of the universe and 
the neutron-to-proton ratio is no longer able to track its equilibrium value: 
it "freezes-out" . 

After this freeze-out, the neutron-to-proton ratio can be approximated 
by 

(I) =^"^'"^^' (^) 

where Tf is the freeze-out temperature, which is determined by setting the 
equality between the expansion and the weak rates. Note that, since the 
weak rates depend on Am, also Tf depends implicitly on it. 

But while reactions (^-(^ are too slow to track the rate of expansion, 
nuclear reactions are still fast enough to keep light elements in NSE and, 
thus, their abundances are still very small. It is not until the third and last 
stage that nuclear production effectively begins. 

• Third Stage: light-element synthesis (0.6 MeV > T > 0.05 MeV; 
3 sec < t < 6 min) 

Shortly before this time, electron-positron pairs have finally annihilated, 
transferring their entropy to the photons and increasing their temperature 
relative to that of the neutrinos by a factor (11/4)^''^. So, at this last stage 
the only relativistic species are the neutrinos and the photons. 

The neutron-to-proton ratio does not remain really constant after the 
freeze-out but it continues to decrease during this stage from its freeze-out 
value because of neutron decay and the effect of the strong neutron-proton 



^The numerical value of the Riemann zeta function is C(3) ~ 1.202. 



reactions. Actually, it will only be at T ~ 0.05 MeV, when practically all 
available neutrons are bound into nuclei, that the neutron-to-proton ratio 
becomes constant. 

But the most important feature of this period is that, while in the previous 
stage only weak processes were relevant, at this epoch strong and electromag- 
netic interactions get important and nucleosynthesis begins. 

The evolution of light-element abundances is dominated by the competi- 
tion between the nuclear reaction rates and the expansion rate. Although the 
densities of the light "fuels" for the reactions involved in the process are now 
significant, they eventually are not big enough to keep up with the demand 
for the NSE of heavier elements. Moreover, Coulomb-barrier suppression 
becomes gradually more important. Both effects result in the freeze-out of 
nuclear reactions and the consequent series of departures of light nuclei from 
their NSE states. 

Let us describe the process in detail, following ||12[. The most abundant 
nucleus, ^He, is to an excellent approximation produced only through the 
mass-3 nuclei, ^He and tritium (t), and it is only through these nuclei and 
the reactions ^He(n, 7)'^He and t{p, 7)'^He that it is allowed to be in NSE. At 
temperatures greater than 0.6 MeV, the ^He and t abundances are sufficiently 
large to allow ^He to track its equilibrium value. But at T ~ 0.6 MeV, 
the NSE curves of the mass-3 nuclei cross with that of ^He, with the ^He 
abundance rising faster than those for ^He and t. At about this temperature, 
also, the reactions that maintain ^He in equilibrium become too slow, ^He is 
forced to leave its NSE curve and follows instead the corresponding ^He and 
t NSE tracks. 

This goes on until T ~ 0.2 MeV. At this temperature ^He and t also en- 
counter a bottleneck, a "minor" deuterium {d) bottleneck: reactions d{n, 'y)t 
and d{p, 'j)t, which keep d in NSE with the mass-3 nuclei slow down and now 
^He, ^He and t follow the d NSE curve. 

At T ~ 0.08 MeV, the ^He(n,p)t reaction freezes-out and the mass-3 
nuclei lose their NSE too. 

Finally, all the elements encounter the "major" deuterium bottleneck at 
T ~ 0.06 MeV and after that the nuclear species evolve in quasi-statistical 
equilibrium. 

Some ^Li and some ^Be are synthesized, but due to the existence of energy 
gaps among stable nuclei at mass numbers A = 5 and A = 8 and significant 
Coulomb-barrier suppression at this time, the production of nuclei beyond 
y4 = 8 is inhibited. 



3 The a-dependence of element abundances 

We are now able to discuss which are the most important a-dependent mag- 
nitudes that affect primordial abundances. 

During the first two stages of SBBN, the only a-dependent parameters 
are the weak reaction rates @-(i|), which in turn determine the the freeze- 
out temperature Tp (see Eq. (^)), the neutron-proton mass difference Am 
(see Eq. (||)) and the binding energies Ba (see Eq. (P)). 

Since the abundances of composite nuclei during these stages are by them- 
selves negligible, the effects of a small variation in their binding energies are 
negligible too. 

But this is not the case with the weak reactions (^-(^ and with the 
neutron-proton mass difference. Since almost all available neutrons after 
freeze-out are finally bound into "^He, in a good approximation its abundance 
is given by 

n - 2-A!^ (8) 

From this equation and Eq. (|^) it follows that the main parameters that 
fix the final value of I4 are precisely Am and Tf. Radiative and Coulomb 
corrected expressions for the weak reactions (|l)-(@) have been calculated in 



13 1 and it was shown that their corrections to primordial abundances are 
very small.^ The effect on Tf of a variation in a is then given only through 



the variation of Am (see [|T3|, Q). That means that Kj depends just on a 
single a-dependent parameter, i.e.. Am. 

This fact makes the constraint on a estimated from the ^He primordial 
abundance very model-dependent: a particular model for the a-dependence 
of Am, is needed to estimate the change in Y4 due to a change in a. One 
may also note that changes in the other gauge coupling constants, such as 
the strong coupling a^, may induce changes that are even more important, 
but more difficult to compute than those of a changing a |jT5|, |16 . 

As an example of the a-dependence of Am, we have used the treatment 



of J. Gasser and H. Leutwyler |T^, from which we deduce 

Am ~ 2.05 - 0.76 ■ (1 + Aa/a) MeV (9) 



The most important correction is that for He, of about 1%. 



This method is a phenomenological way to evolve the mass of the nu- 
cleoli as determined by the change of quark masses due to their strong and 
electromagnetic energy and the electromagnetic binding. It is a qualitatively 
reasonable method but has no fundamental QCD backing. 

Other abundances are not so sensitive to Am as I4: because of its large 
binding energy, ^He acts as a "sink" during primordial nucleosynthesis, its 
abundance is less sensitive to changes in the nuclear reaction rates than the 
other abundances, and more sensitive to variations in the parameters that 
fix the number of neutrons relative to the number of protons, i.e., Am and 
Tf. We note, following |T5[, that the dimensionful weak (Fermi) coupling 



constant Gp does not depend on the gauge coupling constant in the standard 
electroweak model. 

This means that we should use the other light-element abundances to 
limit Aa, in addition to the ^He abundance, bypassing the problem of the 
a-dependence of Am. 

Our next step is, then, to analyze the role of a during the third stage. 
The fine structure constant affects the third stage because it enters into the 
expressions of the nuclear reaction rates: as we have already mentioned. 
Coulomb-barrier penetration is a determining factor during nucleosynthesis. 
To implement the changes in a in these expressions we first need to know 
exactly the way by which thermonuclear reaction rates are obtained from 
experimental data. This will be discussed in the next section. 

4 Thermonuclear reaction rates 

The rate of a nonrelativistic nuclear reaction taking place in a nondegenerate 
environment is as usual given as the thermal average of the product of the 
corresponding cross section (y{E) and the relative velocity times the number 
densities of the particles involved: 

Rij = ninj{a\v\) (10) 

Under the assumption that these particles have isotropic Maxwell-Boltzmann 



kinetic energy distributions this average can be written as |18 



a\v\)=[ — \ T-i Ea{E)e-TdE (11) 

\yU7ry JO 



where /z is the reduced mass. 

For charged-particle induced reactions^ the cross section is given by 

a^E) = ^^e-^-''^^) (12) 

where S{E) is the cross section factor or the astrophysical 5*- factor and ri{E) 
is the Sommerfeld parameter: 



where Eq is the Gamow energy and Zi, Zj the electric charges of the coUiding 
nuclei. 

The reason for introducing this factorization for the cross section is the 
fact that, because of the exponential energy dependence of the Coulomb 
barrier penetrability, charged-particle cross sections are extremely difficult to 
measure at low energies. Therefore, it is necessary to extrapolate cr{E) from 
experimental data to lower energies. Since the exponential factor is given by 
solid quantum mechanical principles, only the unknown nuclear physics part 
S{E) has to be fitted, and it is generally a slowly varying function of energy. 

Let us first consider the non-resonant terms of the reactions 1-12 (see 
tables and |]). Because of its slow variation with energy, we can expand 
S{E) as a Taylor series: 



S{E) = S{0) + ^ 



E + - 



E=o 2 dE^ 



E^ (14) 

E=0 



When inserted into Eq. ([TT|) this gives 
where 

OO _ 1 



E=Vi 



iV^(0= / Z/*e-^e-^^ Uy (16) 

JO 



^We are not concerned with neutron-induced reactions since they do not involve any 
Coulomb barrier penetrability factor in their cross sections and, therefore, they are not 
very sensitive to a pSl]. 



and where 



^ = 27raZa,^. (17) 

Introduce further k = ^^/4. At low temperatures compared with the 
Coulomb threshold, these "astrophysical" integrals are given by |^ 

iVo(0 = 2^Kh-'-'' (l + l.-i) (18) 

iV.(0-2/|.V-(l + |.-i) (19) 

iV,(0=2y|je-^'^^(l + |.-i) (20) 

From these expressions it is straightforward to obtain the a-dependence of 
the non-resonant contributions to the reaction rates we are treating. We have 
not considered the a-dependence of the reduced mass /i and of S{E) and its 
derivatives about the zero energy since this dependence obeys a polinomial 



and, finally 



law on a [|18| and it can be safely neglected. 

For reactions 3, 5, 7 and 8 the non-resonant terms are predicted to have 
a combination of polynomial and decreasing exponential terms |]l8l ; so, ex- 
pressions like 

S{E) = 5(0)g(E)e-^^ (21) 

where q{E) is, again, a slowly varying function of the energy, must be in- 
cluded into the ^'-factor. For most of the cases we are considering, a very 
good approximation results when q{E) is assumed to be a constant, i.e., 
q{E) = g(0) = 1. Then an additional term of the form 

W\v\),, = (—] ' T~lTf,S{0)No{^p) (22) 

where T^ = T/(l + (3T) and ^fs = 2naZiZj{fi/2y^'^T^ ' , appears in the 
expressions of the corresponding reaction rates. 

However, for reaction 3, /? is assumed to be of the form 




/3 = 7o + 7(^l T (23) 

and then the effective temperature T^ is also an a-dependent parameter: 

10 



Tp = — ? (24) 



Continuum and/or narrow and broad resonant terms appear in the ex- 
pressions for the rates of reactions 4-6, 11 and 12. These terms are generally 
of the form 



E 



{a\v\)res = 9{T)e'T (25) 

where E is about the continuum threshold energy or the resonance energy, 
respectively, and g{T) is a function of the temperature. But since the elec- 
tromagnetic contribution to these energies is very small compared with the 
strong contribution, we can safely neglect their a-dependence. 

Finally, reactions 4, 5, 11 and 12 have in addition a cut-off factor for the 
non-resonant terms of the form 

fco = e-^^of (26) 



where Tco is a cut-off temperature. From [|18| it follows that this temperature 
is roughly proportional to a~^. 

In tables ^ and ^ we have explicitly shown the relevant dependence on a 
for all the charged-particle induced reactions having an impact on primordial 
abundance changes. For the relative change in the value of a, we define the 
quantity 5: 

5 Constraints on the variation of the fine struc- 
ture constant 

Using the theoretical framework we have described we proceed to see what is 
the effect of varying a on the relative abundances. We do use the newest ver- 



sion of the SBBN code [|T^ to implement these variations. We keep all other 
coupling constants fixed and assume no other "exotic" effects are present, 
like strong primordial magnetic fields. We discuss these effects below. We 
have performed several tests to see that the code is working correctly. We 
use a neutron lifetime of r^ = 886.7 sec, which is the currently accepted value 
according to the Particle Data Group p|. 

11 



In Fig. |I| we show the results for the final abundances of the observation- 
ally interesting elements ^He, ^He, d and ^Li as a function of 77, for SBBN 
and for a 5 % increase or decrease of the value of a during nucleosynthesis. 
In Figs. ^|-0 the results for the abundances are shown separately on an ex- 
panded scale. As can be seen, the fractional change in the ^He abundance is 
quite insensitive to the value of r/, as are the other abundances with the no- 
table exception of ^Li. This is due to the two competing mechanisms for ^Li 
production, i.e., for 77 < 3 x 10~^°, ^Li is produced by '^He(t, 7)^Li, a process 
in which the Coulomb-barrier is not as significant (but where the change in 
Am caused by a change in a is) as in the reaction ^He(^He, 7)^Be, which, 
followed by decay of ^Be to ^Li, is the dominant reaction that synthesizes ^Li 
for r7> 3 X 10"^°. 

Actually, the of 77 coming from a comparison between the predictions of 

standard nucleosynthesis and observational data is presently quite uncertain, 

as is exemplified by the differing recent analyses of the problem in SBBN 

21, 22|. In |21], a combined analysis gives t] ~ 1.8 ■ 10^^'^ (driven largely by 



accepting a high value of the deuterium abundance), whereas in |2^ a low 
deuterium abundance is taken as the preferred observational result leading 
to a value 77 ~ 5 ■ 10~^°. In Fig. ^ the fractional changes in the light-element 
abundances are shown versus the assumed fractional change in a for both 
these values of tj. 

In the case of the low value of 77 ~ 1.8 ■ 10"^'', the ^Li abundance does not 
vary as strongly with a changing a, and since this abundance is known less 
accurately than that of ^He it is the latter which has to be used to bound 
possible variations of a. However, the observational status ^He is still a 
matter of debate, and it does not seem likely that the systematic errors are 



larger than usually quoted. For instance, the global fit referred to in |2l|] is 



F4 = 0.238 ± 0.002 ± 0.005 (28) 



IS 



whereas the corresponding value in a recent reanalysis by Izotov et al. [^ 
given as 

F4 = 0.244 ± 0.002 (29) 

Using the second upper limit and the first lower limit for the ^He abun- 
dance, as is difficult to avoid before better data become available, we see 
from Fig. |^ (a) that a reduction of the value of a larger than 5 % is not ex- 
cluded, if ^He is considered alone, and the other abundances are only used to 
generously bound t] to be in the interval 10~^° — 10~^. Interestingly, one may 
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note that on the contrary an increase of the electromagnetic strength by 2 % 
would relieve the pressure these low-r^ models feel from the upward revision 
of the ^He abundance of p3[, while still falling in the allowed range for ^Li 



discussed below. Of course, we have to remember the fact that the change in 
'^He is caused by the model-dependent variation of the neutron-proton mass 
difference with a. 

For the high value of ?] ~ 5 ■ 10~^°, we see that a decrease of a by 2 % 
is within our allowed range of the ^He abundance (the change in deuterium 
caused by such a change in a does not affect very much the r] value assigned to 
a given deuterium measurement). Since for the high ri solution the variation 
of ^Li with a is rapid, we have to consider the observational data for lithium. 
Unfortunately, the situation is less than clear here as well. It seems that the 
data are well described by a "plateau" value versus metallicity, given by ||2^ 

Y-j = (1.73 ± 0.30) ■ 10"^° (30) 

The fact that there is little dispersion around this value could indicate 
that it represents the primordial abundance. However, lithium could be 
destroyed by stellar processes by an unknown amount, which has led some 



workers in the field [^] to consider the possibility that the true primordial 
value is up to a factor of 2 larger. 

It is interesting to note from Fig. |1] that the ^Li abundance for this value 
of 7] can be brought down to 2 • 10^^° if the strength of a was larger by 
around 3 % than the standard value (thus making the Coulomb suppression 
stronger). For the same change of a, the "^He abundance is, however, violating 
our (model-dependent) upper bound by a considerable amount. Thus, we 
conclude that the high-rj solution, which already has problems explaining ^Li 
(if ( ppp represents the primordial value) can not be improved by changing a, 
in contrast to the low-r/ solution which may benefit from an increase of a by 
a couple of percent during nucleosynthesis. 

6 Conclusions 

The successful first-order comparison between Big Bang nucleosynthesis pre- 
dictions and observed abundances implies that the values of the fundamental 
coupling constants at that epoch can not have been too much different from 
the present ones. When it comes to more precise statements, the situation 
is complicated by the presently partially conflicting measurements. (Similar 

13 



conclusions have recently been obtained concerning possible limits on the 
effective number of neutrino species during nucleosynthesis |^.) We argue 
that it is impossible to quote a better bound than 2 % on the deviation from 
the present value of a at the time of nucleosynthesis, and note that such 
a variation may even be implied, if the deuterium observations leading to 
the low-r/ solution are confirmed with increased significance, and the trend 
continues of an increasing value for ^He from observations. 

A new ingredient in our treatment is the discovery of the large sensitivity, 
especially at high values of 77, of the ^Li abundance on the electromagnetic 
strength, due to exponential Coulomb barriers. If the experimental and 
theoretical situation is improved for this isotope, it holds the promise of 
giving the most stringent bound on variation of a, with the added virtue 
that it is less sensitive to model assumptions than, e.g., "^He. The problem 
that the high-r^ (low deuterium) solution may have in explaining the low 
plateau value of ^Li can in principle be solved by an increase in a by a few 
percent, but only if the true primordial ^He abundance is well above 0.25 (or 
if Am depends on a in a different way than we have assumed) . 

It is of interest to notice that other exotic effects such as the presence of 
strong primordial magnetic fields mainly affect the ^He abundance ||2^ and 



not Li or the other elements. Also in this respect the signal given by ^Li is 
clearer. 

The limit obtained here is about the same as the one potentially given by 
the recombination time test []TU|, but probably more significant since it is at 



a higher redshift. However, at the time when the microwave background was 
emitted, the universe was totally dominated by gravity and electromagnetic 
interactions from which it can be argued that it offers a less model-dependent 

test. 
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Number 


Reaction 


Reaction rate 


1 


'^H{p,'yfHe 


1 ^ 2 i 

2.650 X 103(1 + 5)3 T„ ^exp(-3.720(l + 5)aT^ '')x 

2 1 2 2 

x(1.000 + 1.120 X 10^1(1 + <5)"3Tg3 + 1.990(1 + 5)^Tg3 + I.56OT9 + 
+1.620 X 10-1(1 + 5)3Tq^ + 3.240 X 10-^(1 + <5)iT|) 


2 


■iH{p,-,)^He 


2.200 X 10'*(l + 5)^T9"^exp(-3.869(l + 5)§Tg"^)x 

x(1.000 + 1.080 X 10-1(1 + <5)-iTg* + 1.680(1 + (5)§r| + 1.260r9 + 

+5.510 X 10-1(1 + 5)3t| + 1.060(1+ 5) ir|) 


3 


^Li(p,jyBe 


6.690 X 10^(1 + 5)iT,lTg^ exp{-8.413{l + S)iT~J) 
Tga =T9/(1.000- 9.690 x lO-^Tg + 2.000(1 + 5)iT|(l. 000 - 9.690 x lO-^Tg)-!) 




4 


<^Li{p,afHe 


3.730 X 10i"(l + 5)iTg"^ea;p(-8.413(l + <5)iTg"^ - (1.818 x 10-iT9)2(l + 5)2)x 

x(1.000 + 5.000 X 10-2(1 + 5)-iTg3 - 6.100 x 10-^(1 + 5)iT^ - 

-2.100 X IO-2T9+ 6.000 X 10-'^(l + <5)3Tg^ + 5.000 x 10-'^(1 + 5)iT9^)+ 

+1.330 X 10"'Tg"*exp(-1.776 x 10iTg-i)+ 
+1.290 X 10''Tq"iexp(-2.182 x lO^T-i) 


5 


'^Li(p,a)^He 


1.096 X 10*^(1 + 5)4rg"^exp(-8.472(l + 5)§Tg"*)- 

-4.830 X 10»(l + 5)3Tg|rg"^exp(-8.472(l + 5)iTg^^)+ 

+1.060 X 10i''Tg"^exp(-3.044 x 10iTg-i)+ 

+1.560 X 10^(1 + 6)iT^^ exp{{-8A72{l + 5)iTg^) - {5.896 x 10-1X9)2(1 + 5)2)x 

x(1.000 + 4.900 X 10-2(1 + 5)-iTg^ - 2.498(1 + 5) irg3 + 

+8.600 X lO-iTg +3.518(1 +<5) 3 Tg^ + 3.080(1 + 5) iT93)+ 

+1.550 X 10<*Tg~2exp(-4.478Tg-i) 
[Tgt =T9/(1.000 + 0.759T9)] 


6 


^H(a,'yfLi 


3.010 X 10^(l + 5)3T,'^^exp(~7A23(l + 5)3T,^^)x 

x(1.000 + 5.600 X 10-2(1 + <5)-iTg^ - 4.850(1 + 5)1 Tg^ +8.85OT9- 

-5.850 X 10-1(1 + 5)3Tg^ - 5.840 X 10-1(1 + 5)^7;*)+ 

+8.550 X 10iTg"2ea;p(-8.228Tg-i) 



Table 1: Reaction rates (in units of cm^s~^mole~^) for the most relevant a- 
dependent reactions involved in SBBN (reactions 1-6). Tg is the temperature 
in units of 10^ Kelvin. 
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Number 


Reaction 


Reaction rate 


7 


'iH{a,^yLi 


3.032 X Vf'{l + S)TiT^^exp(-?,mii(l + 5)iT^^)x 

x(1.000 + 5.160 X 10-2(1 + 6)-iT^ + 2.290 x 10-^(1 + 5)iT| + 

+8.280 X lO-'^Tg - 3.280 X 10-*(1 + <5)lTg^ -3.010 x IQ-'^il + 5)iT^)+ 

+5.109 X 105(l + 5)3T|;T~^exp(-8.068(l + <5)iTg;*) 
[T9c = T9/(1.000 + 0. 13879 )] 


8 


'^*He{a,-yfBe 


4.817 X Vi^{l + 5)^T^^exp(-lA9& x 10^(1 + <5)iT9"^)x 
x(1.000 + 3.250 X 10-2(1 + 5)-iTg3 - 1.040 x IQ-'^ {1 + 5)'^ T^ ~ 
-2.370 X lO-^Tg -8.110 x V)-^'(l + 5)^T^ - 4.690 x IQ-^' (1 + 5)iTh+ 

1 5 3 2 1 

+5.938 X 10S(l + <5)3Tg6^Tg 2exp(-1.286 x 10^(1 + <5) a Tg^^ ) 
[T9d = T9/(1.000 + 0.107T9)] 


9 


'^H{d,nfHe 


3.950 x 10«(l + 5)^T9"^ezp(-4.259(l + 5)§Tg"^)x 

x(1.000 + 9.800 X 10-2(1 + 5)-i7;f +7.650 x 10-^(1 + <5)iTg3 + 

+5.250 X 10-ir9+ 9.610 x 10-^(1 + <5) 3 Tg^ +1.670 x 10-^(1 + 5)1 Tg^) 


10 


'^H{d,pfH 


4.170 x 10»(l + 5)^Tg"^exp(-4.258{l + 5)lTg"^)x 

x(1.000 + 9.800 X 10-2(1 + <5)-iTg3 + 5.180 x 10-^(1 + <5)iTg^ + 

+3.550 X lO-iTg -1.000 x 10-^(1 + <5)lTg3 -1.800 x 10-^(1 + <5)§Tg3) 


11 


^H{d,nfHe 


1.063 x 10"(l + 5)^Tg"^exp(-4.559(l + <5)iT'g"^ - (1.326 x 10iT9)2(l + 5)2)x 
x(1.000 + 9.200 X 10-2(1 + 5)-i7;f -3.750 x 10-^(1 + <5)§Tg^- 
-2.420 X IO-IT9+ 3.382 x 10^(1 + <5)3Tg^ +5.542 x 10^(1 + <5)iTg*)+ 
+8.047 X 10*^Tg"3g^p(_o.486Tg-i) 


12 


^He{d,p)'^He 


5.021 X 10"'(l + 5)^Tg"^exp(-7.144(l + <5)§Tg"^ - (3.704r9)2(l + 5)2)x 
x(1.000 + 5.800 X 10-2(1 + 5)-iT'g^ +6.030 x 10-^(1 + <5)iTg^ + 
+2.450 X IO-IT9+ 6.970(1 +<5) 3 Tg^ + 7.190(1 + 5)iT;f) + 
+5.212 X 10**Tg"^exp(-1.762Tg-^) 



Table 2: Reaction rates (in units of cm^s^^mole"^) for the most relevant a- 
dependent reactions involved in SBBN (reactions 7-12). Tg is the temperature 
in units of 10^ Kelvin. 
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Figure 1: The abundances of the hght elements D, '^He, ^He and ^Li as a 
function of r], the baryon-to-photon ratio. Curves are shown for the standard 
value of the fine structure constant a, and for a variation of ± 5 % of the 
standard value. The results for the individual abundances are shown on an 
enlarged vertical scale in Figs. 
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Figure 2: The abundance of (a) deuterium and (b) "^He as a function of ?7, 
the baryon-to-photon ratio. Curves are shown for the standard value of the 
fine structure constant a, and for a variation of ± 5 % of the standard value. 
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Figure 3: The abundance of (a) ^He and (b) '''Li as a function of t], the 
baryon-to-photon ratio. Curves are shown for the standard value of the fine 
structure constant a, and for a variation of ± 5 % of the standard value. 
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Figure 4: Fractional variation of the light-element abundances Y vs fractional 
variation of the fine structure constant a, for two values of the baryon-to- 
photon ratio 7]. Notice that the curves for ^He for the two rj values nearly 
overlap. The dramatic decrease of the ^Li abundance for a high rj value is 
due to the strong Coulomb barrier in its production. 
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